Local and nonlocal entanglement for quasiparticle pairs induced by Andreev reflection 
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We investigate local and nonlocal entanglement of particle pairs induced by direct and crossed 
Andreev reflections at the interfaces between a superconductor and two normal conductors. It is 
shown theoretically that both local and nonlocal entanglement can be quantified by concurrence and 
detected from the violation of a Bell inequality of spin current correlators, which are determined only 
by normal reflection and Andreev reflection eigenvalues. There exists a one-to-one correspondence 
between the concurrence and the maximal Bell-CHSH parameter in the tunneling limit. 
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Quantum entanglement is a peculiar nonlocal corre- 
lation between distant quantum mechanical systems 1 . 
It plays a key role in the emerging tasks of quantum 
information, such as quantum cryptography 2 , quantum 
dense coding 3 ", quantum teleportations, quantum error 
correction^, and quantum computation^. Entanglement 
is also critical for enabling a quantum computer to per- 
form some computational tasks exponentially faster than 
any classical onei This kind of nonclassical puzzling 
correlation has been demonstrated experimentally for 
photons^ and atoms^. However, it remains a challeng- 
ing task for experimentalists to demonstrate electronic 
entanglement in solids. 

There is growing interest in controllable entanglement 
between electrons in solids, especially in mesoscopic sys- 
tems. A variety of schemes have been proposed, focusing 
on the creation, manipulation, and/or detection of spin 
entanglement of electron pairs by using either quantum 
dotsifi or superconductorsiii^. A mesoscopic normal- 
conductor-superconductor(NS) device has been designed 
for the creation and detection of orbital entanglement be- 
tween electron pairsi 3 -. All of the proposals with hybrid 
NS structures only allow for ideal spin-independent tun- 
neling at the interface between the normal conductor and 
the superconductor. The degree of entanglement is still 
not, to the best of our knowledge, quantified for electron 
pairs emitted from a superconductor. A novel idea for 
electron-hole entanglement produced in tunneling events 
in normal conductors has been recently put forward by 
Beenakker et alii*. Concurrence^ and the Bell-CHSH 
parameters^ for entangled electron-hole pairs have also 
been analyzed based on the scattering matrix theory. 

In this work, we study local and nonlocal entangle- 
ment between pairs of quasiparticles (electrons and holes), 
which are created in Andreev reflection processes in hy- 
brid normal-conductor-superconductor systems. As an 
example, we consider a structure with two normal con- 
ductors connected to a common superconductor via tun- 



nel barriers^, as depicted in Fig. 1. The separation d 
between these two normal conductors should be compa- 
rable to the superconducting coherence length £, in order 
to guarantee the occurrence of the crossed Andreev reflec- 
tion process^. One can expect entangled pairs of cither 
electrons or holes depending on the sign of voltage ap- 
plied to one of the two normal conductors. If a negative 
voltage =FeV is applied to the left conductor, an incident 
spin-up(down) hole (electron) in the left conductor can 
be Andreev reflected as a spin-down(up) electron(hole) 
either (i) at the left conductor, or, (ii) at the right con- 
ductor. The former process is termed as direct Andreev 
reflection, and the latter is crossed Andreev reflectioi 



Since the absence of a hole (electron) in the filled Fermi 
sea is equivalent to the creation of an electron(hole), the 
above process can be viewed as emission into the conduc- 
tors L and R of two electrons(holes) with opposite spins 
from the superconductor. 

We first consider a simple case in which the trans- 
parencies T L / R of the tunnel barriers for spin-up and 
spin-down particles are the same. In the tunnel- 
ing limit (Iz,^ <C 1), the outgoing state will be a 
superposition of the vacuum state |0), the local en- 
tangled state (|T {-E)) L || {E)) L - || {-E)) L || 
(£)}i)/v2, and the nonlocal entangled state (|| 
( -E)) L || (E)) R - II {-E)) L |t (E)) R ) /V2 with weight 
y/l-ri-T L T R X(d), T L ,y/T L T R X(d), where X(d) cx 
e -2d/7rj oiggcj-ikgs distance suppression of breaking-up of 
a pair in s-wave superconductors^. 

Now we consider a general case allowing for an arbi- 
trary value of barrier transparency and spin-flipping at 
the NS interface. Two normal conductors are labeled by 
L and R for brevity. A negative voltage — eV is applied to 
the normal conductor L to create entangled pairs of elec- 
trons. Then the input state describing a stream of holes 
at energies < E < eV injecting to the superconductor 
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FIG. 1: A schematic picture of direct and crossing Andreev 
reflection in the standard semiconductor model: two nor- 
mal conductors (L and R) are connected to a superconduc- 
tor(middle) via tunnel barriers, with a negative voltage — eV 
being applied at the conductor L. The shaded region is the 
filled electron state. A hole(open circle) with spin down(up) 
below the superconducting chemical potential incident on the 
superconductor from L will be either converted as an elec- 
tron(filled circle) with spin up(down) above the chemical po- 
tential at L or at R. The process can be viewed as emission 
of an electron pair from the superconductor. The inset is the 
schematics of a hybrid structure with a superconductor and 
two normal conductors with separation d. 



can be written as 
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(1) 



0<E<eV 



where the ground state \G) is the filled Fermi sea up 
to the superconducting chemical potential in the normal 
conductors. a^ (E) is the creation operator exciting an 
incident spin-<r(=|, I) holc(a = h) or electron(a = e) 
state with energy i5(which is counted from the supercon- 
ducting chemical potential) in normal conductor X = 
L, R, and similarly for operators b of an outgoing state. 

Since in this work we want to discuss the entangle- 
ment arising only from Andreev reflection, a more general 
discussion in the presence of cotunneling events will be 
given elsewhere. The experimental setup should satisfy 
the condition eV <C A (the superconductor energy gap) 
where the cotunneling rate between the two normal con- 
ductors is negligibly small. Taking into consideration the 
crossed Andreev reflection, a scattering matrix relating 
the output state to the input state will be 
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noting the matrix transpose. The 2x2 submatrices S XY 
in Eq. J2J) is spanned in the spin space, and the elements 
of the submatrices are S^y. JiiJ2 . They describe normal 



reflection, direct Andreev reflection and crossed Andreev 
reflection processes. The scattering matrix can be ob- 
tained from the scattering matrix approaches by match- 
ing the wavefunctions at the NS interface. In general, 
the scattering matrix has a pretty complicated relation to 
the elements of the scattering matrices of the barriers in 
the presence of crossed Andreev reflection. Fortunately, 
it takes a simple form in the tunneling limit, as we show 
later. 

Following Beenakker et alii^, we rewrite the input state 
in the following form 
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In Eq. @, <Jy is the Pauli matrix 
fromEq. (2) 



(3) 



One finds 



y=iVs%+b'*sfi.+b'*s% 
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The input state evolves unitarily into the output state 
|*out) = U\V in ) = Hi Ua\W\G) in the Schrodinger pic- 
ture. However, in the Heisenberg picture, state vectors 
remain unchanged, and operators transform as ai «-* bi — 
U^ciiU — ^2jSijdj. The operation Ua\u^ appearing in 
the evolution from the input state to the output state can 
be considered as a time-reversal evolution as compared to 
the original transformation between a^s and 6j's. Such 
an evolution is equivalent to exchanging the role of inputs 
and outputs, i.e., U =>■ U\ S =4* S' and a. L =>■ bi. Then 
we have Ua[U^ =$-Wb\U = J2j Sjibj, which means that 
a direct substitution of b\ 's for a\ 's in the input state via 
the related scattering matrix relation gives the correct 
output state. Therefore, we obtain the output state by 
inserting Eq. (@J into Eq. J3J) 
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This expression reveals distinct physical processes. The 
first term represents normal reflection with the incoming 
holes reflecting backward at the NS interface, the sec- 
ond and third terms describe high-order Andreev reflec- 
tion processes with the incoming holes transforming into 
electrons at the left and/or right conductors, and the last 
term gives direct and crossed Andreev reflection, which 
constitutes the desired local and nonlocal entangled pairs 
of particles. 

Defining some auxiliary states as 
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we rewrite the output state in the following form 

|* OUt > - i{A\^LR) + Ax\lpx) + II (Si h L°yS hhT 
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With the help of unitary conditions of the scattering 
matrix the state normalization condition yields : 
Ax = (TrT X LTl L ) l ' 2 A = (TrTTt) 1 /2 ) with j XL = 

qeh QhhT y _ qeh _ qehT 
D XL a V D LL "y> 1 — LL y RL "V 

To manifest the entanglement more clearly, we trans- 
form from the electron-hole picture to an all-electron pic- 
ture with the particle-hole transformation, c e L \(—E) = 
b h La (E), c\\ Ra {E) = b L \ Ra {E). The new vacuum state 

is defined by |0) = Y[ n<E<eV . a a h L \{E)\G), and the re- 
normalized output state becomes to leading order in the 
Andreev reflection matrix S eh 
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where Cx = A 2 X . The output state is a superposition of 
the new defined vacuum state |0), and local and nonlocal 
entangled states \4>l), \iPr)- The local entangled state 
consists of electron-electron pairs at the same conductors 
and nonlocal entangled states of electron-electron pairs at 
two different conductors. Notice that the entanglement 
states here are a product of two-particle entangled states 
at different energies. 

In the tunneling limit T L / R y^ -C 1, all the elements 
of the Andreev reflection matrices (S XL )ij ^ !• Ex- 
panding the product in the output state iJBJ to first order 
in (S'x L )ij and using the fermion commutation relations 
repeatedly—, the output state after rearranging opera- 
tors becomes 
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Such a state is a sum of superposition states of the 
newly defined vacuum state, and local and nonlocal en- 
tangled states at different energies, which describes a 
wave-packet-like entanglement state. 

For a two-qubit pure state one can quantify the 
degree of entanglement by concurrence^ C = | (i/)\(Ty <8) 
|, where the superscript * denotes complex conju- 
gate. For the local and nonlocal entangled states i/jl's 
and ipRS, one finds 



C^l/r) = 2 



y/ Rl\RliA ll i R fA LL i Rl 
Rl-\A ll / ri + RliA ll / r ^ 



(8) 



which becomes in the tunneling limit (All/r^/^ <C f, 
since T L/m/l < f ) 
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where Rlt, Rli, A ll / r ^ and A LL / R ^ are the eigenval- 
ues of the normal reflection matrix product S L V and 
the Andreev reflection matrix Sll/r^ll/R' respectively. 
In the absence of multiple Andreev reflections, we find 
simple expressions for the eigenvalues of the normal and 
crossed Andreev reflection matrices: A LL y^ T^yjA, 
A ER yi f» X(d)Ti J y^T R yi/4. Maximal entanglement is 
achieved in the spin-independent tunneling case. We em- 
phasize that, the distance suppression of crossed Andreev 
reflection does not degrade the degree of entanglement. 
Its effect is just to diminish the weight of the nonlocal 
entangled state \i/)r). 

Experimentally, the degree of entanglement can be de- 
tected by violation of a Bell inequality^. In solid state 
systems, the Bell inequality can be constructed in terms 
of the spin correlatorsiii&ii 

B = \F(0 X ,0 Y ) + F(0' X ,O Y ) + F(0 X ,0' Y ) - F(O' X ,0' Y )\ 
< 2. (10) 

In Eq. I|1U[1 . the spin correlators are expressed in terms 
of low-frequency spin current correlators 
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where the low-frequency spin current correlators^ are 
obtained as Af** = J d(t - t>) < SI Xai (t)8I Ya2 (t') >= 
e 3 V\(S x h x Sf Y ) ai „ 2 \ 2 x x/h {x = 1 for X ^ Y and 
x = 2 for X = Y). The spin current correlators pro- 
vide information about the joint probabilities to find 
one spin-<7i particle at X along the direction Ox, and 
the other with spin-02 at Y along the direction 0y^. 
The direction along which the detection is performed 
can be incorporated by placing a unitary rotation ma- 
cos Ox — sin 6x 
sin Ox cos Ox 
flection and the Andreev reflection matrices 1 ^. We find 
after some algebra 
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In order to express the Bell- 



1 

-1 

CHSH parameter B in terms of the eigenvalues Rl<t, 
All /Ra as well, we perform the polar decomposition^ of 
the scattering submatrices with these 2x2 unitary ma- 
trices Ux, U 2 , U 3 , V: S h L h L = U X R L V\ S? L = U 2 A LL V\ 
S L h R = UzAl R V\ where the diagonal matrices are = 

diag^R^.R 1 ^), A LL = diag(A L f l v A 1 /^), A LR = 
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diag(A]/jL, A^ 2 , ). Performing a similar calculation as 
in RefsiS 3 - ancU^ results in the maximal Bell-CHSH pa- 
rameter B for the local and nonlocal entangled spin states 
|V>l) and \tps) 



B(iPl/r) = 2 



1 



4:Rl]RliA ll / r ^A ll / ri i 1/2 
(i?L T ^LL/i?T + RliA ll/ri ) 2 - 



(13) 

In the tunneling limit, the maximal Bell-CHSH parame- 
ter becomes 



BiS> L/R ) = 2 



4A LL / R1 A LL/Rl 



(All/ri 



All/ri) 2 



1/2 



(14) 



From Eq. I|14|) one can readily obtain the following 
expected relation B = 2(1 + C) 1 / 2 as compared with 
Eq. 10. The maximal Bell-CHSH paramter always vi- 
olates the Bell inequality as long as Rli, Rli, A ll / r ^ 
All/r.^ 7^ 0- I n the spin-independent tunneling case, 
B = 2y/2, the Bell inequality is maximally violated. As 
expected, there exists a one-to-one correspondence be- 
tween the concurrence C and the maximal Bell-CHSH 
parameter B in the tunneling limit. We notice that, the 
expressions of the concurrence and the maximal Bell- 
CHSH parameter bear formal similarities compared to 
that in tunneling between normal conductors^. 

The dcphasing effect on the entanglement in our sys- 
tem can be discussed in a way similar to Refs^ ancU 3 -. 
Here we discuss possible experimental realizations. For 
pure superconductors, the coherence length is of order 
£ ~ H/Sp ~ hv F /A ~ e F is typically 10 3 - 10 4 

times A, kp is of order of 10 8 cmT 1 , and £ is typically 
10 3 ~ 10 4 nm? A . Therefore, the exponential suppres- 
sion with increasing distance between two tunneling bar- 



riers poses no severe restriction to experimental setup. 
A pair of tunneling barriers has been already prepared 
with distance of the order 10 nmi£. Since kp£ 3> 1, 
the distance suppression is dominated by a polynomial 
term [kpd)~ 2 if the superconductor considered is three- 
dimensional. One way to reduce or exclude such a power- 
law decay is to replace the three-dimensional supercon- 
ductor by a lower-dimensional one 1 ^. Since the Andreev 
reflection process involves two quasiparticles with one 
particle energy above and the other below the super- 
conducting chemical potential, energy-resolving or spin- 
filtering current noise detectors can be used to measure 
the low-frequency current correlators 2 ^. Alternatively, 
one can also use a pair of spin-resolved edge channels in 
the quantum Hall effect proposed by Beenakker et ali 1 ^ 
for our purpose. 

In summary, we consider a unique accessible structure 
in which both local and nonlocal entangled spin states 
for both electrons and holes can be produced at the same 
time. The concurrence and the maximal Bell-CHSH pa- 
rameter have a simple dependence on the eigenvalues of 
the normal and Andreev reflection matrices, and possess 
a one-to-one correspondence in the tunneling limit. Pos- 
sible experimental realization is also discussed. 
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